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Flavourful hadronic physics 
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We review theoretical approaches to form factors that arise in heavy-meson decays and are hadronic expressions 
of non-perturbative QCD. After motivating their origin in QCD factorisation, we retrace their evolution from 
quark-model calculations to non-perturbative QCD techniques with an emphasis on formulations of truncated 
heavy-light amplitudes based upon Dyson-Schwinger equations. We compare model predictions exemplarily for 
the F^~'^{q'^) transition form factor and discuss new results for the go'Dn coupling in the hadronic D* decay. 



1. Flavour physics and strong phases 

The last two decades have witnessed impor- 
tant advances in flavour physics and in partic- 
ular heavy-meson decays. From the first obser- 
vation of a S meson by the CLEO Collaboration 
in 1981 at the Cornell Electron Storage Ring dj 
(and their ongoing I^-meson research program) 
to the dedicated B-physics facilities at SLAC in 
California and KEK in Japan, much progress has 
been made. Of course, while B physics is the 
main focus of the Collaborations Belle at KEK 
and BaBar at SLAC, and of the CDF experiment 
at Fermilab, considerable efforts have also been 
devoted to studies of D-meson decays, charmo- 
nium and r physics. 

Naturally, the driving force is to confirm the 
electroweak sector of the Standard Model which 
has established itself as the foremost paradigm to 
describe CP violation; in other words, the main 
task experimentalists strive for is the precise area 
of the Cabibbo-Kobayashi-Maskawa (CKM) tri- 
angle and the weak CP violating phase codified 
within its angles. In order to measure the size and 
exact form of this triangle, the angles are derived 
from branching fractions in a nowadays large va- 
riety of decay channels. We focus on non-leptonic 
decays, B M1M2, but our discussion also ap- 
plies to the cleaner semi-leptonic B Mtvi case. 

From a theoretical point of view, heavy mesons 
can be used to test simultaneously all the mani- 



festations of the Standard Model, namely the in- 
terplay between electroweak and strong interac- 
tions. They also provide an excellent playground 
to examine non-perturbative QCD effects already 
much studied in hadronic physics. It is notewor- 
thy to remind that no CP- violating amplitude 
can be generated without strong phases. Sup- 
pose a heavy particle H decays into a mesonic 
final state M = M1M2..., H ^ M, and that the 
Standard Model lagrangian contributes two terms 
(two Feynman diagrams) to this process. Then, 
the decay amplitude and its corresponding CP 
conjugate, written most generally, are 
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The weak coupling is a combination of possibly 
complex CKM matrix elements and Ae^''' denotes 
the strong (hadronic) parts of the transition am- 
plitude, where we emphasise that they too can 
have both a real part, or magnitude, and a phase, 
or absorptive part, due to multiple rescattering 
of the final-state quarks and mesons. These CP- 
related intermediate states must contribute the 
same absorptive part to the two decays, therefore 
the strong phases ^p^ are the same in Eqs. ([T]) and 
([2]). Taking the difference of the absolute squares, 
known as direct CP violation, 

I^P - l^p = 2AiA2 Im(AiA;) sin((^i - ^2) , (3) 

one sees that no such violation, |.4/^| 7^ 1, 
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can occur if the weak couplings contain only real 
phases or the strong phases are the same. Hence, 
in order to extract the weak CKM phases with 
precision from the decay amplitudes, it is crucial 
to evaluate the QCD contributions reliably. 

2. QCD factorisation 

As simple as these mesons appear to be — 
a bound colourless heavy-light qq pair — the dif- 
ference in quark masses and the energetic light 
mesons produced in their decays lead to a dis- 
parate array of energy scales. A central as- 
pect of heavy-meson phenomenology are factori- 
sation theorems which allow for a disentangle- 
ment of short-distance or hard physics, which en- 
compasses electroweak interactions and perturba- 
tive QCD, from long-distance or soft physics, gov- 
erned by non-perturbative hadronic effects. In 
the following, we illustrate the factorisation with 
non-leptonic decays of a heavy meson H. 

In the hamiltonian formulation of heavy-quark 
effective theory (HQET) 2 , in which amplitudes 
are expanded in powers of Aqcd/ti/i and the 
heavy quark is a static particle in the limit to/j 
cxD, the H — > M1M2 decay amplitude is given by 



G 



(4) 



where A. 



p 
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ucts of CKM matrix elements and Gi? is the Fermi 
coupling constant. The dimension-six four-quark 
operators Oi result from integrating out the weak 
gauge bosons in the operator product expan- 
sion and the Wilson coefficients Ci(C) encode per- 
turbative QCD effects above the renormalisation 
point C- 

In what is called "naive" factorisation, the 
hadronic matrix element {MiM2\Oi\H) is ap- 
proximated by the product of two bilinear cur- 
rents, (Mi |g7^(l-75)6|B) ® {M2\q'-^^,{l-l^)q\G) 
-\- (Ml M2), where colour indices have been 
omitted. This factorisation simply expresses the 
matrix element of a local four-quark operator as 
a product of a decay constant and a transition 
form factor. However, as has long been known, 
the saturation by vacuum insertion fails in the 
case of most D decay modes and is largely in- 



sufficient to reproduce the experimentally more 
precise data on i? ^ M1M2 branching fractions. 
In fact, any hard final-state gluon interaction has 
been neglected and soft-gluon exchange is at best 
incorporated into an effective colour parameter or 
form factors. Moreover, the renormalisation scale 
and scheme dependence of Ci (C) are not cancelled 
by those of the matrix element {MiM2\Oi\H){C,) . 

A major improvement over this simple fac- 
torisation Ansatz comes from the systematic re- 
organisation of weak and QCD interactions in 
the HQET. Three distinctive approaches have 
emerged in recent years: QCD factorisation 
(QCDF) % perturbative QCD (pQCD) 4 and 
soft-coUinear effective theory (SCET) |5J. We 
here focus on the form factors that emerge in 
these factorisation approaches and solely remark 
that QCD corrections beyond naive factorisation 
entail, in the limit rrih 3> Aqcd, the B M1M2 
decay amplitude can be schematically written as 



(M1M2IO.IB) - (Mi|ji|B)(M2b2|0) 

l + ^r„a^+C'(AQCD/mfc) 



(5) 



where ji and j2 are the bilinear currents. This 
has been shown explicitly to leading order in 
cts [3|6J and including the one-loop correction 
(a^) to the tree-diagram scattering between the 
emitted meson and the one containing the spec- 
tator quark [7]. 

The factorisation theorem derived using SCET 
agrees with QCDF if perturbation theory is ap- 
plied at the hard rri^ and hard-coUinear rrifcA 
scales, with A typically of the order of 100 MeV. 
It is evident from Eq. ([5]) that higher orders in as 
break the factorisation but these corrections can 
be systematically supplemented; the analogy with 
perturbative factorisation for exclusive processes 
in QCD at large-momentum transfer is not acci- 
dental [5]. Further contributions that break the 
factorisation, formally suppressed in AqcD/^h 
yet not irrelevant, are weak annihilation decay 
amplitudes and final-state interactions between 
daughter hadrons [Qj. Neglecting power correc- 
tions in as and taking the limit mb — > 00, the 
naive factorisation is recovered. 
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3. Separating scales: the softer the harder 

While factorisation theorems elaborated with 
SCET provide the means to systematically inte- 
grate out energy scales in the perturbative do- 
main, yielding approximations valid in the heavy- 
quark limit for a given decay in terms of products 
of hard and soft matrix elements, a reliable eval- 
uation of the latter is notoriously difficult. In 
fact, it is the soft physics of the bound states 
that renders the task hard, as it implies non- 
perturbative QCD contributions. Full ab initio 
calculations are currently out of reach and for the 
time being one is left with modelling the heavy- 
to-light amplitudes with as much input from non- 
perturbative QCD as possible. Just how much 
soft physics is included depends on the theoreti- 
cal Ansatz and techniques employed. 

A variety of theoretical approaches have been 
applied to this problem, recent amongst which are 
analyses using light-front and relativistic quark 
models, light-cone sum rules (LCSR) and lattice- 
QCD simulations. In Section [4] we briefly sum- 
marise these approaches while Section [5] deals 
in more detail with studies of heavy-to-light 
form factors within the framework of the Dyson- 
Schwinger equation (DSE). We refer to a recent 
review |10) for a summary of transition form fac- 
tor data from lattice-regularised QCD and just 
note that current results are obtained at large 
squared momentum transfer, i.e., ~ 16 GeV^ 
for i? ^ TT transitions. Hence, values at low 
must necessarily be extrapolated by means of a 
(pole-dominance) parametrisation |ll|23j . 

4. Hadronic transition form factors 

Quark models: relativistic quark models f l2ll3l 
114115116117118] have in common that their only 
degrees of freedom are constituent quarks whose 
masses are parameters of the hamiltonian. The 
hadronisation of the two valence quarks is de- 
scribed by vertex wave functions or phenomeno- 
logical Bethe-Salpeter amplitudes (BSA). The 
approaches in [13ll4ll5j represent heavy-to-light 
transition amplitudes by triangle diagrams, a 3- 
point function between two meson BSA and the 
weak coupling, which yields the transition ampli- 
tude {M {p2)\qT jh\H (pi)) and reads generally. 



A{pi,P2) ^tl'CD J J^^^M\k;-'P2)Sq{k+p2) 

X Ti{p,,p2)SQ{k+pi)TH{k;pi)Sg,{k) , (6) 

where S{k) are quark propagators, Q ~ c,b] 
q = q' = u,d, s; M = S, P, V, A and the index ^ 
indicates a possible vector structure in the final- 
state BSA. F/ is the interaction vertex whose 
Lorentz structure depends on the operator Oi in 
the HQET and Th is the heavy meson BSA. The 
trace is over Dirac and colour indices. 

An analogous approach represents the ampli- 
tude in Eq. ^ by relativistic double-dispersion 
integrals over the initial- and final-mass vari- 
ables pI and P2, where the integration kernel 
arises from the double discontinuity of the tri- 
angle diagram (putting internal quark propaga- 
tors on-shell via the Landau-Cutkosky rule). The 
meson- vertex functions are given by one-covariant 
BSA [16117] . Other quark models [12] represent 
i? — > M form factors by overlap integrals of me- 
son wave functions, obtained from confining po- 
tential models, and appropriate weak interaction 
vertices. Similar quark model calculations were 
performed on the light cone [18] . 

All these approaches model soft contributions 
with vertex functions, while the propagation of 
the constituent quark, S{k) — — mq)~^ , is 
scale independent and does not describe confine- 
ment and dynamical chiral symmetry breaking 
(DCSB). As noted in Refs. [171191201 . this can 
lead to considerable model dependance at larger 
momentum transfer. 

Light- cone sum rules: In a LCSR the operator- 
product expansion of a given correlation function 
is combined with hadronic dispersion relations. 
The quark-hadron duality is invoked: the corre- 
lator function is calculated twice, as a hadronic 
object and with subhadronic degrees of freedom. 
After separation in HQET of the heavy meson's 
static part, Ph = p + q = rnuVh + k, where Vh 
is the four-velocity and k is the residual momen- 
tum, and likewise redefinition of the momentum 
transfer q = mhVh + q^p + q = k, one obtains 
the heavy-limit correlation function, 

H^(p,g) = II^^(p,g) +0(l/m,.); (7) 
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g2 [GeV2] 


Ref. [12] 


Ref. [m 


Ref. [H] 


Ref. [20] 


Ref. [2? 


Ref. El 


Ref. [21] 





0.217 


0.29 


0.247 


0.24 


0.25 ±0.05 


0.26 ±0.03 


0.27 ±0.02 ±0.07 


10 


0.41 


0.54 


0.49 


0.53 


0.54 


0.51 




15 


0.67 


0.84 


0.97 


0.91 


0.90 


0.81 




20 


1.40 


1.56 


> 10 


1.75 


1.83 


1.58 





Table 1 

Numerical comparison for the transition form factor, Ff^'^{q^); the values of Ref. [2U) are calculated, whereas 
for Refs. [12I16I18I22I23) the value F^^'" (0) and the corresponding extrapolation in these references are employed. 
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whereas a hadronic correlator can be written 
)\JM\M{p)){M{p)\Jh\HiPH)) 



, (8) 



where J/i(0) and J/i„(0) are heavy-light currents 
and Jm{x) the interpolating current for a pseu- 
doscalar or vector meson. 

In Eq. ([5]) only the light-meson contribution 
is represented but higher and continuum states 
can also be taken into account. In Eqs. ([T]) and 
([5]) , the usual role of the correlation functions has 
been reversed |21I22| : the correlation function is 
taken between the vacuum and the on-shcU B- 
state vector using its light-cone distribution am- 
plitude (DA) expansion and the pion is interpo- 
lated with the light-quark (axialvector) current 
Jm(x). The i?- meson DAs are universal non- 
perturbative objects introduced within HQET. 

In Ref. [33] J however, the correlation function 
is taken between the vacuum and the light-meson 
state, whereas the B meson is interpolated by 
the heavy- light quark current J/i(0). As a re- 
sult, the long-distance dynamics in the correla- 
tion function is described by a set of light-meson 
{tt, K, p, K*) DA. In the last step, a Borel trans- 
formations is apphed to both, Eqs. ((Tj) and (O, 
from which one derives a transition form factor 
expressed as a sum rule. The transformation in- 
troduces a scale via the Borel parameter which, 
in turn, is fixed with sum rules for light-meson 
decay constants. 

Besides a systematic uncertainty owing to the 
duality assumption, the main incertitude lies 
within the DAs which encode the relevant non- 
perturbative effects. Only their asymptotic form 
is known exactly from perturbative QCD. As of 



yet, the first two Gegenbauer moments of the DA 
for various light pseudoscalar and vector mesons 
have been obtained from QCD sum rules with 
very large errors, though the first moment is con- 
sistent with lattice calculations [H]. Moreover, 
the transition form factors must be extrapolated 
to space-like momenta. 

For purpose of comparison, we list _B — > tt tran- 
sitions form factors, F+{q'^), for various models in 
Table [TJ As observed therein, there is a 30% vari- 
ation within the quark models [12|16|18] at = 
which increases at larger q^ values. The LCSR 
predictions [22123] agree at = but their re- 
spective slopes for q^ > vary by 12%. The form 
factors obtained with the DSE model Ref. [20] 
are calculated on the entire physical momentum 
domain and the chiral limit is directly accessible. 

5. Flavourful Dyson-Schwinger equations 

The elements entering the amplitude in Eq. ([6]) 
can be motivated by the solutions of DSEs ap- 
plied to QCD. A general review of the DSEs can 
be found in Refs. [25 26^ and their applications 
to heavy-light transition form factors have been 
investigated in [20127128] . 

Dressed quark propagator: The mesons are 
bound states of a quark and antiquark pair, where 
for a given quark flavour their dressing is de- 
scribed by the DSE (in Euclidean metric), 

S-\p)^Z2{ij-p + m'"^^) + ^ip') , (9) 

with the dressed quark self energy, 

E(p2) = Zi52^V''(p-fc)^7^5(fc)r«(fc,p),(10) 

where represents a Poincarc invariant regu- 
larisation of the integral with the regularisation 
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mass scale A. The current quark bare mass m^™ 
receives corrections from the self energy S(p^) in 
which the integral is over the dressed gluon prop- 
agator, D^^(k), the dressed quark-gluon vertex, 
r°(fc,p), and A" are the usual SU(3) colour ma- 
trices. The solution to the gap equation Q reads 

S{p) = -ij ■ p av{p'^) + crsip^) 

= [z-f ■ p Aip') + B{p^)]-\ (11) 

The renormalisation constants for the quark- 
gluon vertex, Zi{(, A^), and quark- wave function, 
Z2(C,A^), depend on the renormalisation point, 
C, the regularisation scale. A, and the gauge pa- 
rameter, whereas the mass function M(p^) = 
B{p'^)/A{p'^) is independent of Since QCD is 
asymptotically free, it is useful to impose at large 
spacelike the renormalisation condition, 

S-\p)\p2=^2=tj-p + miC^) , (12) 

where to(C^) is the renormalised running quark 
mass, so that for 3> ^qqd quantitative match- 
ing with pQCD results can be made. 

Infrared dressing of light quarks has profound 
consequences for hadron phenomenology : the 
quark-wave function renormalisation, Z{p'^) = 
1/A{p'^), is suppressed whereas the dressed quark- 
mass function, A/(p^) = B{p'^)/A{p'^), is en- 
hanced in the infrared which expresses dynam- 
ical chiral symmetry breaking (DCSB) and is 
crucial to the emergence of a constituent quark 
mass scale. Both, numerical solutions of the 
quark DSE and simulations of lattice-regularised 
QCD [30], predict this behaviour of M(p^) and 
pointwise agreement between DSE and lattice re- 
sults has been explored in Refs. [31 32^. Studies 
that do not implement light-quark dressing run 
into artefacts caused by rather large light-quark 
masses |17I19| to emulate confinement. This is 
because unphysical thresholds in transition am- 
plitudes can only be overcome with the prescrip- 
tion that rriH < rUq^ -\- rriq^ , which poses problems 
for a description of light vector mesons {p,K*), 
heavy flavoured vector mesons {D*,B*) and for 
P-wave and excited charmonium states. 

Whereas the impact of gluon dressing is strik- 
ing for light quarks, its effect on the heavy quarks 
is barely notable. This can be appreciated, for in- 
stance, in Fig. 1 of Rcf. 28J: for light quarks. 



mass can be generated from nothing, i.e., the 
Higgs mechanism is irrelevant to their acquiring 
of a constituent-like mass. 

Bethe-Salpeter amplitudes: The BSA can be 
determined reliably by solving the Bethe-Salpeter 
equation (BSE) in a truncation scheme consistent 
with that employed in the gap equation ([9]) . Con- 
sider the inhomogeneous BSE for the axialvector 
vertex Fg^ in which pseudoscalar and axialvector 
mesons appear as poles: 

rlf^ik; p) = - a^jy^^'^^ - 

X Sf{q+)Til{q-P)Sg{q^)^Tl{q^,k^) (13) 

+ jy-P{k~q)^^^Sf{q+)^kilp{k,q;P). 

In Eq. p^ . P is the total meson momentum, 
q± ^ q± P/2, k± = k± P/2, A^'^^ is a 4-point 
Schwinger function entirely defined via the quark 
self energy f33^ and /, g denote the flavour indices 
of a light-light or heavy-light bound state. The 
solutions of the vertex F^^ must satisfy the axial- 
vector Ward-Takahashi identity, 

P^ Vilik- P) = Sj\k+)l^r, + ll,S-\k^) 

-zK(C)+mg(C)]F^^(fc;P), (14) 

where F^^ solves the pseudoscalar analogue to 
Eq. (jl3p . A systematic, symmetry-preserving 
truncation of the DSE and BSE is given by the 
Rainbow ladder [34] which is their leading-order 
term with the dressed quark-gluon vertex, F^, re- 
placed by 7;^. It can be shown that A^^^ = in 
this approximation. 

The above Ward-Takahashi identity possesses 
another remarkable property; the set of quark- 
level Goldberger-Treiman relations that follow 
from it reveal the full structure of the light-meson 
BSA [35^. In particular, it enables one to relate 
the leading covariant, e^{k\ P), of the light pseu- 
doscalar BSA with the scalar part, B{p^), of the 
dressed-quark propagator (|lip in the chiral limit. 
This motivates an effective parametrisation of the 
light mesons (M — tt, K), 

Fm (fc; P) = ^75 £M (fc') = 175 Bm {k'')/.fM , (15) 
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/m = Im I V^, which has been capitahsed on in 
transition form factor calculations 20,27,28] . 

Simultaneous solutions of the quarks 's 
DSE and the heavy meson's BSA with 
renormalisation-group improved ladder trunca- 
tions, obtained for the kaon |35| . prove to be 
difhcult. The truncations do not yield the Dirac 
equation when one of the quark masses is large. 
A recent attempt to calculate BSA for D and B 
mesons reproduces well the respective masses 
but underestimates experimental leptonic decay 
constants by 30 — 50%. With a consistent deriva- 
tion of the heavy meson BSA pending, simple 
one-covariant forms for TH{k;P) are currently 
employed in Eq. ([5]), which reproduce leptonic 
decay constants in a simultaneous calculation. 

6. Hadronic decays 

The decay D* — > Dtt can be used to extract 
the strong coupling g between heavy vector and 
pseudoscalar mesons to a low-momentum pion in 
the heavy meson chiral lagrangian |37| . One con- 
siders the matrix element, 

{D{p2)7T{q)\D*{pi,X)) ^go'Dn ex-q, (16) 

where the couphng, go'D-n ~ 17.9 ± 0.3 ± 1.9, is 
experimentally known 38J and related to g. Sim- 
ilarly, one may also extract g from the unphysical 
decay B* — > Bn in the chiral limit and where 
mb/AqcD corrections are better controlled. 

The coupling gwD-n is related to a heavy-to- 
heavy transition form factor via the LSZ reduc- 
tion of the pion and the use of PCAC, 7r(x) — 

Ai^i{x) / {fTr'm'^) , which leads to: 

{D{p,)niq)\D*ip,)) = q^ ^'^J-f^ 

X J d'xe^^-^Dip2)\A^{x)\D*{p,)). (17) 

Hence, the matrix element in Eq. (fT6|) has been 
reduced to the Fourier transform of a transition 
matrix element between the D* and D mesons 
in the chiral limit with the axial QCD current 
Af^{x) — (f'^fi'^^q'' ■ Results from this reduction 
procedure have been obtained on the lattice [321 
I40| and most recently in a simulation with rif = 
2 [H] which yields go-Dix = 20 ± 2. 



This form factor can also be calculated 
straightforwardly without reduction of the pion 
employing Eq. ^ with the dressed quark prop- 
agators in Eq. (fTT|) and substituting the pion's 
BA (lisp for F/. In a reassessment and improve- 
ment of a calculation of go'D-n within a Dyson- 
Schwinger model [28], we obtain go'Dn = 21 [42] 
in agreement with the lattice result [H] and about 
16% larger than the experimental value. 

7. Conclusive remarks 

We have stressed the importance of hadronic 
effects in decays of heavy-flavoured mesons and 
portrayed the various theoretical Ansdtze for the 
heavy-to-light transition form factors. In short, 
the main obstacle to their precise calculation, 
which veraciously reproduces the infrared fea- 
tures of QCD, are the uncertainties of the light- 
cone DA in the case of LCSR and the lack of 
model-independent wave functions in relativis- 
tic quark model calculations. We have argued 
that the running quark mass of the DSE quark 
propagators is crucial to include confinement and 
DCSB effects in the transition amplitudes; an un- 
finished task are consistent solutions of the BSE 
for the D and B mesons within the DSE formal- 
ism, which will reduce model dependence. 
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